Bright Ultraviolet Regions and Star Formation Characteristics in 

Nearby Dwarf Galaxies^ 



Nicholas W. Melena^ 
Lowell Observatory, I4OO West Mars Hill Road, Flagstaff, Arizona 86001 USA 

nmelenaOlowell . edu 
Bruce G. Elmegreen 

IBM T. J. Watson Research Center, 1101 Kitchawan Road, Yorktown Heights, New York 

10598 USA 

bgeSus . ibm. com 

Deidre A. Hunter, Lea Zernow^ 
Lowell Observatory, I4OO West Mars Hill Road, Flagstaff, Arizona 86001 USA 
dahSlowell . edu , IzernowShmc . edu 

ABSTRACT 

We compare star formation in the inner and outer disks of 11 dwarf Irregular 
galaxies (dim) within 3.6 Mpc. The regions are identified on GALEX near-UV 
images, and modeled with UV, optical, and near-IR colors to determine masses 
and ages. A few galaxies have made 10^ — 10^ Mq complexes in a starburst phase, 
while others have not formed clusters in the last 50 Myrs. The maximum region 
mass correlates with the number of regions as expected from the size-of-sample 
effect. We find no radial gradients in region masses and ages, even beyond the 
realm of Ha emission, although there is an exponential decrease in the luminosity 
density and number density of the regions with radius. Ha is apparently lack- 
ing in the outer parts only because nebular emission around massive stars is too 
faint to see. The outermost regions for the 5 galaxies with HI data formed at 
average gas surface densities of 1.9-5.9 Mq pc~^. These densities are at the low 
end of commonly- considered thresholds for star formation and imply either that 
local gas densities are higher before star formation begins or sub-threshold star 
formation is possible. The first case could be explained by supernovae triggering 
and other local processes, while the second case could be explained by gravita- 
tional instabilities with angular momentum loss in growing condensations. The 
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distribution of regions on a log(Mass) — log(age) plot is examined. The distribu- 
tion is usually uniform along log(age) for equal intervals of log(Mass) and this 
implies a region count that varies as 1/age. This variation results from either 
an individual region mass that varies as 1/ age or a region disruption probability 
that varies as 1/age. A correlation between fading-corrected surface brightness 
and age suggests the former. The implied loss of mass is from fading of region 
envelopes below the surface brightness limit. 

1. Introduction 

Dwarf irregular galaxies are among the most numerous galaxies known and as such are 
a fundamental component of the universe. They arc the most common type of galaxy within 
the Local Group, and are thought to be the building blocks of larger galaxies. Despite their 
numbers and proximity they are still among the least understood of objects. We don't know 
how these galaxies form stars or how they have evolved to become the galaxies we see today. 
As such, it is essential that we work towards understanding the forces at work within them, 
and the processes which have shaped them. 

Understanding the large and small-scale star formation processes in dwarf galaxies has 
long posed a problem to the current models. In the classic Toomre instabihty (Toomre 
1964) there is a specific column density (Ec) needed for a ring-like perturbation to grow 
exponentially. This model seems to work well in the inner regions of most spirals (Kennicutt 
1989; Martin & Kennicutt 2001). Dwarf galaxies are so bereft of gas, however, that in many 
cases Sc is never reached, even in the centers where the gas column density is highest (unless 
the velocity dispersion is very low - de Blok & Walter 2006). Still, we see star formation 
not only in these centers, but even in the outer reaches of the disk where the gas density 
to critical density ratio (Ep/Ec~ 1/Q) is considerably lower than unity. Similarly, the outer 
disks of spirals and some entire spiral disks are subcritical (Martin & Kennicutt 2001). One 
possibihty is that localized processes dominate in the absence of large-scale instabilities (e.g., 
Elmegreen & Hunter 2006). Another possibility is that stars and gas together promote the 
gravitational instabilities that lead to star formation (e.g., Yang et al. 2007). A third 
possibility is that the gas is unstable even in subthreshold environments, but the instability 



^ Based on observations made with the NASA Galaxy Evolution Explorer ( GALEX). G ALEX is operated 
for NASA by the Cahfornia Institute of Technology under NASA contract NAS5-98034. 

^Current address: The University of Arizona, 1401 E University Boulevard, Tucson, Arizona 85721 USA 

^Current address: Harvey Mudd College, 301 Piatt Boulevard, Claremont, California 91711 USA 



-3- 



growth time is long - comparable to the time for angular momentum to leave a growing 
perturbation through viscous and magnetic forces (Elmegreen 1987). In all three cases, 
dwarf galaxies offer insight into star formation without the influence of spiral density waves. 

In this paper we compare the properties of individual star-forming regions in outer and 
inner disks of a sample of dwarf galaxies. The young regions are identified on GALEX 
(GALaxy Evolutionary eXplorcr; Martin et al. 2005) near-UV (NUV) images. We use far- 

UV (FUV) and NUV, UBV, and some JHK colors to model the masses and ages of the 
regions. Our investigation includes approximately 2000 separate star-forming regions. With 
these data we compare mass, age, and several other key properties of star-forming regions 
in inner and outer disks. 

Star forming regions in galaxies are also useful for studies of cluster evolution, partic- 
ularly regarding fading from stellar evolution and mass loss, and cluster disruption through 
cloud collisions and other processes. An important diagram for these studies is a plot of 
log M versus log t for cluster mass M and age t. Such a plot shows a uniform distribution of 
points over a range of log t for each interval of log M above the magnitude detection thresh- 
old for the survey (Hunter et al. 2003; Rafelski & Zaritsky 2005; Fall, Chandar & Whitmore 
2005; Bastian et al. 2005a; Whitmore, Chandar, & Fall 2007). The plot also shows a cluster 
count decreasing as 1/M over a range of log M for each interval of log t. The implication of 
the former trend is that cluster count decreases in equal intervals of age as dN/dt ocl/t 
for all observable masses (Fall et al. 2005). The implication of the latter is that A^ decreases 
in equal intervals of mass as dN / dM cx 1/ for all observable ages (Elmegreen & Efremov 
1997; Zhang & Fall 1999; Hunter et al. 2003; de Grijs & Anders 2006). Here we study these 
diagrams for our sample of galaxies, finding the same trends, and also study in more detail 
the distribution of region surface brightness and radius with age. 

In §2 we discuss the images used in our survey. §3 describes the means by which we 
calculated the photometry for our star-forming regions, as well as a discussion of our age and 
mass modeling process. One of the principle goals of this survey was to examine radial trends 
within the dwarf galaxies, and §4 discusses these radial trends. In §5 we present a discussion 
of some of the major results of our survey: §5.1 considers the surface brightness evolution of 
regions in four best-case galaxies, showing a trend with age that is steeper than from stellar 
evolution fading alone. This suggests a specific pattern of mass loss in each region that results 
from surface brightness limitations. It is a mass loss that is independent of stellar evolution, 
evaporation, and region disruption. Size-of-sample effects regarding surface brightness are 
also discussed. §5.2 reviews the log mass — log age plots with a comment about the implied 
time dependence of region mass, §5.3 suggests that outer disk Ha is missing because of 
faintness, not variations in the initial mass function (see also Hunter, Elmegreen & Ludka 
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2009), §5.4 comments on variations in the star formation rate, §5.5 examines outer disk HI 
column densities and suggests mechanisms for sub-threshold star formation, and §5.6 notes 
the size-of-sample effects in our observations. Lastly, §6 provides a general overview and 
summary of what we have learned about the star formation processes in dwarf galaxies. 
Also included is an appendix in which the properties of the individual galaxies within our 
survey are discussed. 

2. The Images 

The galaxies in this investigation are a sub-sample of a large multi-wavelength survey 
of 136 relatively normal nearby galaxies without spiral arms (94 dim, 24 Blue Compact 
Dwarfs, and 18 Sm galaxies; Hunter & Elmegreen 2004, 2006). That survey is representative 
of the range in galactic parameters exhibited by dwarfs, and is aimed at understanding star 
formation processes in dwarf galaxies. The UBVJHK and Ha data obtained as part of that 
survey are described by Hunter and Elmegreen. 

To extend this survey into the UV, we obtained images from the GALEX archives for 
44 of our 136 dwarf galaxies. These data are part of a study of azimuthally-averaged UV 
properties in Hunter et al. (2009). For the present study, we selected 11 dim galaxies from 
this GALEX sample, all within 3.6 Mpc of the Sun, to measure the properties of individual 
star-forming regions. These data are combined with comparable measurements on our UBV 
and Ha images and on JHK images of one of the galaxies. The galaxies that are part of 
the present study are listed in Table [1] along with several basic properties. 6*74 Li^X exposure 
times and tile names are given in Table [H 

GALEX images were taken simultaneously in two channels: FUV, with a bandpass of 
1350-1750 A, an effective wavelength of 1516 A, and a resolution of 4.0" and NUV with a 
bandpass of 1750-2800 A, an effective wavelength of 2267 A, and a resolution of 5.6". The 
images were processed through the GALEX pipeline, and we retrieved final intensity maps 
with a 1.5" pixel scale from the archives. The GALEX field of view is a circle with 1.2° 
diameter; we extracted a portion around our target galaxies. 

We edited foreground stars and background galaxies, and removed sky from the GALEX 
NUV and FUV images. In some cases the sky was a constant determined from regions 
around the galaxy, but in other cases the sky was determined from a low-order, two- 
dimensional fit to the surroundings. We geometrically transformed the UV images to match 
the orientation and scale of our UBVJHK images. In that way we could identify regions on 
the NUV images and then easily measure fluxes on all of the UBVJHK, FUV, and NUV 
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images using the same extraction apertures. 

3. Photometry and Modeling 

We identified individual star-forming regions on the NUV images, measured their fluxes 
on all of our broad-band images, and subtracted the background galaxy light. We corrected 
these fluxes to the standard photometric systems. We then apphed Bruzual & Chariot (2003) 
stellar population models to our colors and luminosities to determine the ages and masses of 
the star-forming regions. These ages and masses were studied with respect to galactocentric 
radius and other variables. 



3.1. Region Selection 

We marked with polygons all of the individual knots of star formation seen in the NUV 
images. These regions are defined as spots or clumps that stand out against the diffuse 
background disk emission. The polygons were then applied to all the other image passbands 
that we have for each galaxy. 

The regions were identified by eye. One source of uncertainty comes from the occasional 
difficulty of distinguishing foreground stars from star-forming regions in the galaxies at the 
resolution of the NUV images, especially for the more distant galaxies. Although foreground 
stars are preferentially red, we did not wish to bias the regions to only the youngest by 
selecting only blue objects. Instead, wc limited the sample to galaxies within 3.6 Mpc where 
the resolution of the NUV image is about 100 pc. Thus, most young regions are resolved. 

3.2. Background Subtraction 

The star-forming regions sit on a diffuse disk of galactic UV emission. We needed 
to determine the underlying emission for each region in order to measure the photometric 
properties of the star-forming sites themselves. Where we could, we subtracted the mode of 
an annulus starting 5 pixels out from the longest radius of the defining polygon and extending 
another 5 pixels in radius. 

However, because of the clustered nature of these sources, this was out of the question 
for the majority of them. There were often one or several other regions within bounds of a 
reasonable annulus centered on the target region. We specifically wanted to avoid subtracting 



- 6 - 



one region from another, so in these cases we found several background spots close to, but 
separated from, the star-forming region and used the median in these regions to characterize 
the underlying disk. The separated background apertures identified on NUV images were 
also used as background regions for the other passbands. The background emission was 
defined to be the median pixel value within the background region closest to the star-forming 
region. This median pixel value was subtracted from each pixel in the target star-forming 
region to get the flux from the star-forming region alone. 

To check our determination of the background in the case of the separate background 
regions, we determined the sky using an alternate algorithm (the mode) in one galaxy, DDO 
50, using the same background regions we had used with the original (median) photometry. 
The purpose was to ensure that the background value was not biased by the occasional star 
that would appear in a background region in the optical that was not there in the UV. 
The lack of disagreement or the presence of an offset between background values determined 
using the two methods suggest that our background values in separate regions are robust 
against this particular problem. 



3.3. Photometry 

After the underlying galactic disk was subtracted, we transformed the fluxes to the 
standard GALEX, UBV, and JHI^ systems. The UV data are in AB magnitudes. The 
UBV photometry is on the Landolt (1992) system, and the JHK fluxes were calibrated with 
UKIRT standard stars (see Hawarden et al. [2001] for a discussion of near-IR photometric 
systems). We calculated IsAnuvi My, and as well as the colors FUV — NUV, B — V, 
U — B, J — H, and H — K. The uncertainties in the magnitudes and colors reflect Poisson 
statistics in both the signal and in the sky. 

We corrected the photometry for extinction using the same total reddening E(i? — V)t 
as was used in our earlier study of the UBV JHK images (Hunter & Elmegreen 2006). The 
total reddening E(i? — V)t is foreground E(i? — V)f from Burstein & Heiles (1984) given 
in Tabled] plus a constant 0.05 mag for reddening internal to the dwarf galaxy. A constant 
internal reddening of this order for the stars is consistent with measurements of the Balmer 
decrement in HII regions in a sample of 39 dim galaxies (Hunter & Hoffman 1999). There 
the average reddening in HII regions is 0.1, and we have taken half this to represent the stars 
outside of HII regions. We combined the E(i? — V)t with the reddening law of Cardelli et al. 
(1989) to produce the extinction in each filter. For the GALEX?Ateis, Apuv = 8.24:E{B-V)t 



^JHK data were only available for one galaxy in our data set, NGC 2366. 
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and Anuv = 7.39E(i? — V)t- The photometry for all of the regions are given in Table [3l 
which is available in its entirety on-line, and the JHK photometry for the regions in NGC 
2366 are given in Table HI 

There are, in total, 1920 regions for all 11 galaxies. There are FUV — NUV colors 
for 1913 of them; 1792 of the 1920 have B — V colors and 43 have J — H. We removed 
regions from our sample with various problems with the photometry, including photometric 
errors greater than 1 magnitude, contamination by foreground stars, unphysical colors, and 
incomplete photometry. We were left with 1623 star-forming regions with good photometry. 

One reason for the incompleteness of optical and near-IR colors is that star-forming 
regions bright in the NUV are not necessarily bright in the other bands, and can be lost in 
the galaxy background. Regions with a few thousand counts in the UV might only have a few 
hundred counts in V . In addition, the lower the flux that a region has, the more susceptible 
it is to uncertainties in sky subtraction. For any of our regions with highly uncertain colors, 
there was a good probability that our models would be unable to find a realistic fit. That is 
why we removed all of the star-forming regions with any photometric errors greater than 1 
magnitude. 

Another factor working against the smaller regions is that they are often too faint to be 
seen at greater distances. While our survey is complete at a level of —10 or —9 mag in My, 
below an My of —7 mag it is incomplete. For the closer galaxies, those less than 1 Mpc away, 
such as NGC 6822 and IC 1613, this was not a problem, but at the longer distances of > 3 
Mpc this effect became more apparent. This problem is illustrated in Figure [H where the 
distributions of My for all of the star forming regions in NGC 6822 and DDO 50 are shown. 
DDO 50 at 3.4 Mpc has a higher proportion of brighter regions, probably from blending 
effects, and NGC 6822 at 0.5 Mpc peaks at much fainter levels. 

To explore the consequences of distance on our results, we have taken WLM, which 
is nominally at a distance of 1 Mpc, and smoothed the images with a Gaussian kernel 
to the resolution we would expect at a distance of 3.4 Mpc, the distance of DDO 50. We 
then analyzed the smoothed images, including identifying UV-bright regions and background 
regions, performing the photometry, and running models. The results are shown in Figure 
[2] where we compare the region ages and masses as a function of radius from the original 1 
Mpc images and from the pseudo-3.4 Mpc images. At 3.4 Mpc, we see fewer regions with 
more mass and shallower radial gradients, as expected, but general trends remain. Age-mass 
ratios, for example, are roughly the same. From a similar comparison of the other figures that 
we use in our discussion, we believe that the results for the more distant galaxies included 
here, particularly DDO 50 at 3.4 Mpc, are as valid as those for the closer systems, especially 
since our emphasis in this study is on radial trends within galaxies, rather than differences 
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between galaxies. 

3.4. Mass and Age Modeling 

In order to compare star forming regions, we determine the ages and masses by finding 
the best model fit to the colors for each region. The model colors are from the Bruzual & 
Chariot (2003) model stellar populations with the "Padova 1994" stellar evolutionary tracks 
(Alongi et al. 1993, Bressan et al. 1993, Fagotto et al. 1994a,b). We use a Chabrier (2003) 
stellar initial mass function (IMF). This IMF adopts the Salpeter (1955) IMF for stellar 
masses above 1 Mq . Below 1 Mq the IMF is shallower than Salpeter's and below 0.3 
Mq , it flattens out. Mass-to-light ratios arc of order 1.4-1.8 times smaller with a Chabrier 
IMF compared to a Salpeter IMF. Wc used the Z = 0.004 or Z = 0.008 mctallicity models, 
whichever was closer to the metallicity expected from nebular oxygen abundances (Hunter & 
Hoffman 1999) or, for DDO 210 and LGS3, estimated from the mass-metallicity relationship. 
Thus, we used Z = 0.004 models for all galaxies except NGC 2366 and NGC 6822. We also 
use an exponentially decaying star formation history with a decay time of 1 Gyr for all 
regions. For the short ages we obtain here, this means nearly continuous star formation over 
the age of the region. 

Our modeling program loops over age and finds model colors. For each model color we 
subtract the observed color and divide the difference by the uncertainty in the measured 
color. We then sum the square of this difference over the colors. This is a x^. Among 
all possible solutions, we then average with an exp(— 0.5%^) weighting factor to obtain a 
final age and mass. The uncertainties in many of our colors are sufficiently small compared 
to the model spectral energy distribution (SED) fits that the values were unreasonably 
large (often greater than 40). So, we multiplied all of the color uncertainties by a constant 
factor of 5 for all regions in order to get non-zero Gaussian weights. This is effectively the 
same as considering greater than statistical uncertainties, or systematic uncertainties, in the 
observations and models. In a small number of cases the method does not find a good 
solution at all. We considered this to be the case when 5 or fewer of the age models out of the 
164 trial ages in the Bruzual & Chariot tables with ages larger than 2.3 Myr gave < 40. 
We then modified our method to obtain the best possible SED fit by averaging together the 
trial ages and masses from only those age trials that had the lowest rms deviation between 
the observed and modeled colors (i.e., not dividing by the measurement errors as in the 
method). The weights for this average were taken to be Gaussian functions of this rms 
deviation, rather than x^. The uncertainties in the results for these cases were taken from 
the rms values in the averages. The final model age, mass, and rms of the fit are given for 
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each region in Table [31 

As a test of the modehng, we tried our model on some well known clusters: NGC 104 
(47 Tuc), a Galactic globular cluster, and NGC 604, a large young OB association in M33. 
We obtained UBV colors for NGC 104 from Harris (1996) and measured them for NGC 604 
from the Local Group Survey Data (Massey et al. 2006). We measured FUV — NUV for both 
clusters from GALEX archive images. Both sources had their background light subtracted. 
We adopt an E(5 - V) of 0.04 for NGC 104 from Salaris et al. (2007) and Harris (1996), and 
use Z = 0.004 models, the closest model to the metalhcity of the cluster ([Fe/H]= —0.76, 
Harris 1996). We find an age of 12 ± 1 Gyrs, and a mass of (8.9 ± 0.9) x W^Mq for NGC 
104. Our age agrees very well with recent estimates: 10-13 Gyrs (Salaris et al. 2007), 
12.5±0.5 Gyr (Liu & Chaboyer 2000), 10.9±1.4 Gyr (Salaris et al. 2004), and 11.5±0.8 Gyr 
(VandenBerg 2000). Similarly, our mass agrees well with the mass of 7 x IO^Mq determined 
by Meylan (1988) from dynamics, surface photometry, and a two power-law IMF of NGC 
104. 

For NGC 604, based on the number of observed O stars and a Salpeter stellar IMF, 
the mass was previously determined to be 8 x IO^Mq, with an age of 3-5 Myrs (Hunter et 
al. 1996). For the models, we use an E{B — V) of 0.05, as adopted by Hunter et al., and 
a Salpeter stellar IMF. The oxygen abundance of NGC 604 is about half solar (Diaz et al. 
1987, Caffau et al. 2008), so we use the Bruzual & Chariot (2003) tables with 0.4 times solar 
metallicity. Then, including the GALEX FUV — NUV colors, the models give a mass of 
(1.8 ± 0.1) X IO^Mq and an age of 5 Myrs. These agree with measured values to within a 
factor of 2. The model age and mass increase by factors of 2 and 1.3 ii B — V and U — B 
colors are increased and decreased, respectively, by 0.1 mag to account for bright emission 
lines in the U and V bands. 

Because our survey was based on NUV images, it was not sensitive to red, and therefore 
old, things. It is doubtful that we would see modest star-forming regions much older than 
500 Myrs because of this. There is also some uncertainty in the models due to the rapidly 
changing AGB phases for the age range of 0.3-2 Gyrs (Maraston 2005). While this problem 
is more pronounced in the IR than it is in the optical, between these two problems we feel 
that our models are reasonably accurate up to ~500 Myr of age, and this is the age range 
we emphasize here. 

Similarly, we have a lower age limit of 2.3 Myr. This is the point where the colors begin 
to change with population age. To start earlier than this is to say that we can actually see 
a cluster that is 10^'^ yr old, for example, and in fact, such clusters are probably obscured 
in the UV. There is no color and stellar evolution before this time in the models. 
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One additional issue merits discussion. In modeling star clusters in the Antennae galax- 
ies, Fall et al. (2005) found a sharp-edged gap at ages of 10 Myrs. This is due to the integrated 
light of massive clusters at high metallicity being dominated by red supergiant stars. The 
cluster colors change so abruptly during this phase (8 - 16 Myrs) that the model ages become 
degenerate and avoid ages of 10 Myrs. This is easily understood in the "Conti scenario" of 
massive star evolution in which massive stars spend proportionately less of their He-burning 
lifetimes as RSGs rather than WRs at higher metallicities (Massey & Johnson 1998). On the 
other hand, we benefit from the dependence of massive star evolution on metallicity as this 
effect will be less pronounced in low metallicity systems (compare, for example, Z = 0.02 
with Z = 0.004 in Figures 57 and 59 of Leitherer et al. 1999). Furthermore, for Fall et 
al. this issue introduced errors of less than 0.4 in logt. Thus, any gaps we see in the age 
distributions here are interpreted as lulls in the cluster formation rate. 



4. Observed Radial Trends 

With this survey we now have masses and ages for 1623 star-forming regions in 11 local 
group dims. We also have a rough measure of sizes from the polygon fits. Over the range of 
our survey we see ages from our lower limit of several million years to our completeness limit 
of 500 Myrs, and masses from millions of solar masses down to just dozens of solar masses. 
The individual galaxies are discussed in Appendix \M 

We are particularly interested in how different properties of the star-forming regions 
vary with radius within each galaxy, specifically mass, age, size, and number of regions per 
unit area. To show a few of these basic properties, for each galaxy we plot mass against age, 
mass against galactocentric radius, and age against radius in Figure O The figure shows 
several things: the relative abundances of star-forming regions at various masses and ages, 
the lower mass sensitivity limit for each galaxy as a function of age, the upper age limit 
for consistent detections, and the lower age limit of the models. In the plots of mass or 
age against radius, it is possible to see radial trends that may be present. We also show 
the history of star formation in equal intervals of log-age in Figure HI This plots the total 
number of regions, and the number of regions more massive than the observational limit at 
an age of 500 Myr, and the maximum mass of the regions, all as functions of age. 

We are interested in comparing the properties of star-forming regions in the inner and 
outer disks. We have made this division in three ways. First, we cut at the radius Rna 
beyond which Ha emission is no longer detected. In some of our galaxies we find a significant 
number of star-forming regions past Rho, and we compare the regions inside Rho with those 
outside Rho- These are shown as different colors in Figures [3] and HI Generally there is little 
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difference in the masses and ages for the regions inside and outside Rna- 

Second, we cut by the number of regions we identified, comparing the inner 50% of star 
forming regions to the outer 50%. Galaxies with few UV-selected regions beyond Rna are 
displayed this way in Figure [3] as indicated by the panel labels. We show this also for NGC 
6822, DDO 50, and NGC 2366 in Figure [5l where it is clear for NGC 6822 that there is little 
difference in the masses and ages for these two sectors. This is typically the case. DDO 
50 and NGC 2366 differ, however, in that the masses are systematically lower in the outer 
regions, as shown in Figure [3] also. 

Finally, we cut at the radius where there is a break in the broad-band surface brightness 
profile Rsr- Three galaxies (DDO 75, DDO 216, and M81dwA) show double exponential 
radial profiles: sharp changes in the rate of decline of the exponential disk seen in broad-band 
stellar surface brightness profiles. These breaks are seen in optical and mid-IR (3.6 and 4.5 
/xm) passbands that emphasize all but the youngest stars and do not represent the "edge" 
of the disk, but rather a change in the stellar profile. The breaks do not correspond to a 
particular surface brightness limit, and their relationship to identified star-forming regions 
is inconsistent. The break in DDO 216 occurs well past where we see star-forming regions; 
in M81dwA there were only two identified young regions anyway; and in DDO 75 there are 
a substantial number of regions past the break. For DDO 75, we compare the ensemble of 
regions on either side of the broad-band profile break in Figure [6l showing little difference. 

The last two galaxies in Figure E] are plotted without distinction between inner and 
outer regions: LGS3 and M81dwA have a few UV regions but no Ha regions. 

For galaxies with enough star-forming regions, we look at the luminosity and number of 
regions per unit area as functions of radius. These are shown in Figure [71 The distributions 
are similar to each other because the region masses, luminosities, and sizes do not vary 
much with radius. Figure [8] shows the summed angular filling factors of the UV star- forming 
regions versus radius and Figure [9] shows the region sizes (the radius of a circle that has the 
same area as the polygon that defines the region) measured on the NUV images. Recall 
that the sizes are only estimates from the subjective polygonal outlines of the regions; the 
filling factors use these sizes, and are therefore only approximate as well but for looking for 
general trends they are adequate. The figure indicates that filling factors approach ~ 20% in 
the inner regions. The filling factor profiles are similar to the luminosity and count density 
profiles because the sizes are about constant with radius. Exceptions occur for the few cases 
where region size or mass suddenly change at large radii, as for DDO 50, where the mass 
decreases and the size increases beyond 4 kpc, as compared to the rest of the galaxy. In this 
case the filling factor, which follows the size and the count density, drops more slowly in the 
outer part than the luminosity density, which follows the mass and the count density. This 
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trend seems to be an exception, however. 

The five types of radial profiles — region mass, luminosity density, count density, filling 
factor, and size — consistently indicate that the exponential profile of star formation, and 
presumably the exponential disk itself, is the result of an exponentially decreasing count of 
star-forming regions, each with about the same mass, luminosity, and size. The filling factor 
of star formation decreases too. The interpretation is that the basic star-forming unit has 
constant mass and luminosity throughout a galaxy disk, even beyond the edge of visible HIT 
regions, while the star formation rate is proportional to the number density of these regions, 
which happens to be proportional to the local stellar surface density. The reason for the 
latter proportionality is unknown, but has the effect of reinforcing the existing underlying 
exponential profile with each addition of new stars. 

There is generally a smooth continuation of radial trends beyond Rna, which means 
that the loss of Ha does not correspond to a change in any physical process related to 
star formation. It is more likely an observational artifact, resulting perhaps from an HII 
region surface brightness that is below the detection limit in the outer galaxy. This is to be 
expected in a low pressure environment (Elmegreen & Hunter 2006), because of both the 
diffuse nature of ionization bounded HH regions at very low pressure, and the likely escape 
of Lyman continuum photons from the star-forming regions when the surrounding density 
is low. 



5. Discussion 

5.1. Surface Brightness Evolution 

Surveys of star-forming regions are limited by surface brightness for the detection of 
faint extended objects, rather than by magnitude for point sources. To understand our 
detection limits better and to investigate the evolution of surface brightness, we determined 
the average NUV surface brightness of each region, /inuv, from the NUV magnitude on the 
AB scale interpolated to an arcsec angular size. For resolved regions, this surface brightness 
is independent of galactic distance. It has a lower limit that depends on the exposure time 
of the image and perhaps the de-projected galactocentric radius (because of the bias that 
a selected region must stand out above the surrounding disk). The lower limit should not 
depend on age. 

Figure dU] shows the NUV surface brightness versus galactocentric distance for four 
of our galaxies that have a large sample of regions. The lower limit to surface brightness 
depends on exposure time (Table 2): NGC 6822 has over 3.4 times the exposure time in the 
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NUV than the other galaxies, and a hmit that is ~ 0.7 mag arcsec"^ fainter. There is no 
evident dependence on distance from the galaxy center in /xnuv for NGC 6822, but there is 
a small dependence for the other galaxies in the figure. This reflects a limited bias in the 
selection of regions for measurement, but should not be important in studies involving age 
or mass correlations. 

Figure [11] shows /iNuv versus age for the same four galaxies. Dashed lines indicate the 
fading trend from the evolution of a single stellar population model, for which luminosity 
decreases with time approximately as t~^'^^ ( §5.21) . Star-forming regions approximately follow 
the fading trend down to the surface brightness limit, although the distribution of points 
is actually steeper than this for some galaxies. This steepness suggests either a physical 
expansion over time for each region or a loss of mass over time, beyond the loss from stellar 
evolution (which is included in the fading trend). The oldest regions are all above the fading 
trend because of surface brightness limitations. That is, many old regions are missing from 
our survey, i.e., those below the surface brightness limit on the extrapolated downward trend 
in the figure. We have not lost proportionally as many younger regions by surface brightness 
fading. The dividing age is around 10* yrs for this sample. In general, the dividing age 
between detection and loss from surface brightness limits depends on the exposure time and 
the background noise in the image (in addition to the local galactic surface brightness, as in 

Fig. m- 

There is an interesting size-of-sample effect in Figure [TTl The upper limit to the surface 
brightness decreases with age much more slowly than the lower limit, indicating that most 
of the older regions still observed today were intrinsically brighter than the younger regions 
when they formed. Only the faintest observed old regions, i.e., those at the surface brightness 
limit, are an extrapolation of young regions along the fading trend, and then it is only the 
brightest young regions that can be extrapolated to fit in this way. This increase of initial 
surface brightness with age is a size-of-sample effect, considering that the sample size of 
regions in logarithmic intervals of increasing age is also increasing in proportion to the age. 
The older logarithmic bins include more total regions and therefore a few rare regions of 
intrinsically large initial surface brightness. Extrapolated back to their youth, the brightest 
old regions today would have had /xnuv ~ 21 mag arcsec"^ at an age of 10^ yrs. There are 
no regions in our galaxies as bright as this today, although the embedded bright source in 
NGC 2366 would be of comparable brightness in the absence of extinction. 

Figure [12] shows the region radius versus age for the same four galaxies as in Figure [TT] 
The radius is about constant with age, suggesting very little expansion. Considering that 
surface brightness is from a combination of mass, luminosity-to-mass ratio, and radius, this 
constant radius implies that the steep decrease in surface brightness with age seen in Figure 
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[TT] is from a combination of stellar evolution fading and a loss of stellar mass. This mass 
loss must be in addition to supernovae and stellar winds, i.e., to normal evolutionary effects. 
The implication is that detected clusters have lost mass since they formed. 

Figure [13] shows the fading-corrected absolute magnitude of each region in these four 
galaxies, versus the log of the age, t. Fading correction is done by subtracting (0.69/0.4) logt 
from each absolute magnitude. Recall that the luminosity of a region is proportional to the 
region mass, M, multiplied by the light-to-mass ratio, and this latter quantity varies with 
age approximately as t"*^'^^ ( §5.2p . The luminosity is also proportional to 10~*'-^*^nuv^ gg it 
follows that -0.4Mnuv = const. + log M - 0.69 logt, or 

(-1/0.4) M + const. = Mnuv - (0.69/0.4) logt (1) 

The quantity on the right is shown in Figure [T^ versus logt. The dense clusters of points that 
were sloping downward steeper than the fading line in Figure [TT] are now sloping upward in 
Figure [131 which is the direction of fainter fading-corrected magnitudes over time. A dashed 
line with a slope of +2.5 guides the eye. A general slope f3 on this diagram corresponds to 
a correlation logM = const — 0.4/3 logt from equation [H As (3 ^ 2.5 fits the distribution of 
points in the dense portions of Figure[Tn]for WLM, DDO 50, and perhaps IC 1613, the region 
masses there seem to be decreasing as M oc t~^ over and above stellar evolution effects. We 
do not see this trend in NGC 6822. This type of mass-age relationship will also appear in 
§5.21 in reference to the distribution of points on the log M — log t plots in Figure [3l and to 
the near constant counts of above-threshold regions in the age histograms of Figure |H 

Also in Figure [131 there is a long tail of points at large age with a slight downward slope. 
These are the overly bright old regions also discussed for Figure [Til Their distribution was 
flat on that diagram because they paralleled the constant surface brightness limit, and so 
here they decrease as (—0.69/0.4) logt by construction. We are not interested in them for this 
diagram, but only in the increasing trend among the denser distribution of points, which 
had the puzzling result in Figure [TTl that it was declining faster than possible with pure 
evolutionary fading. 

5.2. Time Dependence of Cluster Mass 

Figures [31 and m showed the number of massive star- forming regions formed as a function 
of time over the past 500 Myrs. At the upper end of the range of ages in our survey, even the 
most massive clusters begin to drop out because of fading. The shift in minimum detectable 
mass goes approximately as Mjade = 982(t/Gyr)'^'^^ Mq, so this minimum mass at 500 Myrs 
is 5 times larger than at 50 Myrs (Hunter et al. 2003). Figure [31 shows this limit by plotting 
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the minimum mass at 500 Myrs as a horizontal dashed hne, and the slope of the mass-age 
relationship as a solid green line. The observed minimum mass for each age is comparable to 
the detection limit. The upper end of the green line in Fig. [3] is positioned at the intersection 
of the horizontal dashed line and the age of 500 Myrs. 

Clusters are increasingly lost from our survey over time because of fading effects as the 
main sequence turn-off climbs down the mass function, and because of cluster disruption, 
stellar evaporation from the cluster, and surface brightness limits. Piskunov et al. (2006) 
find the average open cluster lifetime in the Milky Way to be 322 ± 31 Myrs. Boutloukos & 
Lamers (2003) and Gieles et al. (2006) suggest that the disruption time for a 10"^ M© cluster 
depends on the galactic environment, including the tidal field, molecular cloud collisions, 
and shear rate. Boutloukos and Lamers find that the disruption time in the SMC is very 
long — 8 Gyr, much longer than the time scale they derive for spiral galaxies, and longer than 
the time derived by Piskunov et al. for the Milky Way. Hunter et al. (2003) also suggest 
that clusters in the LMC and SMC have long disruption times. 

The number of clusters per unit age in each mass interval decreases like ~ 1/ age, which 
in a mass-age diagram like Figure [31 means that the density of points is approximately 
constant along horizontal strips (Fall, Chandar, & Whitmore 2005). This constancy is also 
evident in Figure [Has a nearly flat histogram of cluster counts above the threshold mass limit. 
Such a time depletion imphes that clusters are losing mass or getting disrupted constantly. 
In addition, the mass function of clusters in the LMC remains a power law for all observed 
times (Elmegreen & Efremov 1997; Hunter et al. 2003), with no preferred loss of low mass 
clusters in older age groups. In the case where each cluster loses mass like M oc 1/t, which 
is one way to give the observed distribution on the mass-age plane, the lack of a low mass 
turnover in the cluster mass function implies that the lowest mass clusters are lost from the 
survey because of fading limitations, and not disruption. Because the galaxies in the present 
survey are dwarf irregulars with internal tidal forces like those in the SMC and LMC, and 
because brightness limits are higher for the present survey than for the LMC and SMC, we 
expect that fading is the dominant cause for the loss of low mass clusters in the present 
survey too. If physical disruption were the cause instead, then there would be a gap between 
the lowest mass cluster for each age on the mass-age plane, and the fading limit. There is 
no such gap in Figure [31 

In another interpretation of cluster evolution, clusters keep their mass for a time t 
and then get disrupted suddenly, perhaps in a CMC collision. In this case, the disruption 
probability per unit time should decrease with time as 1/t to get a constant number per unit 
log t in fixed intervals of mass, as suggested by the distribution of points on the mass— age 
plane. 
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Our observations cannot distinguish between these two mechanisms of physical cluster 
loss, i.e., between mass loss from each cluster with cluster mass oc 1/t, or sudden cluster loss 
with a disruption rate oc 1/t. In both cases, the loss mechanisms have to be independent 
of cluster mass if the mass function has the same power law slope for all ages. The first 
mechanism cannot be simple evaporation because evaporation gives a constant cluster loss 
rate, dM/dt =constant for a constant cluster density (as determined by tidal effects, for 
example; Baumgardt & Makino 2003). To have M oc 1/t for each cluster, we need an 
evaporation rate dM/dt oc 1/t^ or dM/dt oc M^, which is not the way evaporation works. 
The second mechanism cannot be simply cloud-cluster collisions in the usual model because 
the survival probability in a random ISM scales with time as exp(— t), from Poisson statistics. 

The discussion in §5.11 suggests that some individual region masses decrease with time 
as 1/t. There may be some sudden disruption of regions too, with probability oc 1/t, but 
we cannot see that in our data (we know nothing about when the regions that have been 
disrupted were disrupted). 

A schematic model for mass loss by fading below the surface brightness limit of a survey 
is shown in Figure [TH The two Gaussian curves represent projected profiles I through two 
clusters having the same total mass. One cluster is bright and the other is dim and slightly 
larger in radius. The horizontal red line is the surface brightness limit, which is above the 
outer envelopes of the clusters. The green curves are the running masses measured only 
above the surface brightness limit. They increase from left to right as the cumulative cluster 
masses increase from left to right in the projected images. These running masses are given by 
J^^ 'Kxl'{x)dx for position x and starting position xq on the left in the diagram. The intensity 
in the integral, /' is taken to be zero for the parts of the cluster below the surface brightness 
limit, and the actual cluster intensity / for parts of the cluster above the surface brightness 
limit. The dispersion of the bright cluster is 1 on this scale, and the dispersion of the dim 
cluster is 1.682, designed to make the two clusters have the same size out to the surface 
brightness limit. This satisfies the constraint suggested by Figure [12] that clusters do not 
appear to expand much with age. In fact, if we identify the faint cluster in Figure [14] as an 
old cluster, then it has expanded with age, but at the same time as it has faded, leaving the 
apparent size at the surface brightness limit about constant. The point of Figure [T^ is that 
two clusters with the same total mass can have two different measured masses if only those 
parts of the cluster that are above the surface brightness limit are measured. Cluster fading 
then automatically causes a measured mass loss over and above stellar evolution effects. 
This mass loss is not the same as stellar evaporation because the missing stars are still near 
the cluster and maybe even bound to the cluster. They are just in the faint, unobserved, 
envelope. 
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Several of our galaxies have peaks in the number of star-forming regions as a function of 
time, as well as peaks in the number of massive regions and in the maximum mass of clusters. 
Prom a study of integrated Ha equivalent widths of a large sample of dwarfs, Lee et al. (2007) 
conclude that dwarf galaxies should experience 1-3 starbursts per Gyr. These bursts should 
have a duration of 50-100 Myrs and an amplitude of 6-10 times the star formation rate of the 
quiescent period. McQuinn et al. (2009) have found large-scale bursts of duration 200-400 
Myrs in several nearby dwarfs, with smaller scale "flickering" on timescales of 3-10 Myrs. In 
our sample we see two gaps of duration 60 Myrs and 1300 Myrs over the past 1 Gyr in NGC 
6822. The enhanced cluster formation episode centered on an age of 1 Gyr, for example, 
lasted about 470 Myrs. Thus, the history of cluster formation in NGC 6822 has been roughly 
consistent with the findings of these studies. However, in 6 other systems in our sample, 
cluster formation appears to have been roughly constant until 9-50 Myrs ago, with few or 
none formed since then. Three galaxies show no obvious gaps in cluster formation over the 
past 500 Myr. 



5.3. Comparing Inner and Outer Galactic Disks: The Lack of Ha 

We see no consistent change with radius in the average mass or size of star-forming 
regions for each galaxy. Regions beyond Rua have masses as high as (1.8 ±0.9) x lO^M© for 
the globular cluster in WLM, and (3±1) x IO^Mq beyond Rna in NGC 6822. This lack of a 
mass gradient is consistent with the lack of any changes in the UV radial surface brightness 
profiles at the point where Ha ends (Boissier et al. 2007). Exceptions are DDO 50 and 
NGC 2366, which show the clearest gradients in region mass. DDO 50 and NGC 2366 are 
compared to a more standard case, NGC 6822, in Figure O Mass and age histograms are 
plotted for each galaxy, with the inner and outer half of the regions separated for comparison. 
The median log (mass) for NGC 6822 is only 0.12 larger in the inner half than the outer 
half, while the median log mass is 0.25 larger in the inner parts of DDO 50 and NGC 2366 
than the outer parts. The median ages are slightly younger in NGC 6822 and NGC 2366 
and slightly older in DDO 50 in the outer halves. 

Ha becomes invisible in the outer parts of galaxy disks even though we find normal 
star formation there. For DDO 210, there is no Ha emission even though there is one 
region younger than 10 Myr. There are many reasons why Ha should disappear when star 
formation continues, including a lack of massive stars due to old cluster age or a steep 
IMF (see, for example, Pfiamm-Altenburg & Kroupa 2009), an escape of Lyman continuum 
photons from the vicinities of the 0-type stars, and a undetectable surface brightness for the 
ionization-bounded HII region. The age gradient explanation is difficult to reconcile with 
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our observation of a nearly constant age with radius, and that, along with possible IMF 
variations, are also in conflict with the observation of a normal exponential UV profile. The 
UV exponential implies that the onset of star formation in the outer disk is in the same 
proportion to the surrounding stars as in the inner disk, and in a steady state, this implies 
that the evolution of 0-stars off the main sequence is in the same proportion too. Thus the 
fraction of star-forming regions in each stage of evolution should be independent of radius, 
as should the IMF, and this means that the fraction that produce Ha is constant too. An 
exception might occur if the outer regions had a burst of star formation at just the right 
time in the past to have all of the associated 0-type stars now evolved off the main sequence. 
This is unhkely for all of the galaxies here. 

Escape of Lyman continuum radiation from the vicinity of the 0-star can explain the 
lack of Ha whether or not those photons escape the galaxy. Ha emission disappears when 
the ionizing radiation from a massive star spreads out over a large volume and the emission 
measure becomes undetectable. Considering that the quantity nl^^Re ^ 50 cm~^ pc for 
an 0-type star (Vacca et al. 1996) with Hll region electron density Ue and radius Rg, and 
that the density in the ambient interstellar medium should be 0.1 cm~^ or lower in the 
outer regions of most galaxies, a typical Hll region radius would be Re > 180 pc and the 
corresponding emission measure would be nlRe < 1.8 cm~^ pc, which is undetectable in 
our survey by a factor of 3-40. The hmit of our detection for emission measure is 5.3-74 
pc cm~^ (Youngblood & Hunter 1999), which corresponds to a minimum density of 0.4-1.4 
cm~^ at this ionization rate. Lower densities, particularly off the midplane, would make Rg 
even larger than several hundred pc, and then the Lyman continuum photons could escape 
the galaxy completely. 



5.4. Trends in Recent Cluster Formation Histories 

Our dim galaxies fall into one of three categories describing the overall recent star 
formation histories of the past 500 Myrs: 

Star formation with gaps. These galaxies take two forms. First, NGC 6822 is an example 
of periodic spikes in star formation over a finite duration, and it has made regions as massive 
as ~ 3 X 10^ Mq. The most recent gap in cluster formation occurred around 100 Myrs ago 
and lasted about 60 Myrs. The second form are those galaxies with an apparent cessation 
in cluster formation in recent times. These galaxies include DDO 50, DDO 70, DDO 216, 
and WLM. The recent gap began 9 (DDO 60 and DDO 75) to 25 (WLM) Myrs ago. DDO 
70 and DDO 75 sit in between these two types: They appear to have had a fairly constant 
level of star formation for the last several hundred million years until about 10 Myrs ago. 
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In the past 10 Myrs only a few clusters have formed. In addition, clusters older than about 
200 Myrs are missing except for a few especially massive ones: clusters up to 4 x 10^ Mq in 
mass. 

Constant star formation. These galaxies appear to have had a reasonably steady, 
nonzero rate of star formation over the last 500 Myrs. They make a 5 x 10^ Mq region 
on the order of once every 10 Myrs. They make a 5 x IO^Mq mass cluster about once ev- 
ery 100 Myrs, and some have managed a IO^Mq region at some time in their recent past. 
Galaxies in this category include IC 1613 and NGC 2366. 

Nearly no star formation. These galaxies show up as little smudges barely above the 
background sky on our GALEX images. LGS3 is our archetype for this category, with only 
6 useable regions in the last Gyr, all only a few 100 Mq . We also include DDO 210, as we 
found only 9 useable regions, none of which were more massive than ~ IO^Mq; none of the 
5 from the last 500 Myrs was more massive than 10^ Mq. We also include MSldwA in this 
category; this galaxy has only two measured regions. These three galaxies are also the only 
galaxies in our survey that show no measurable HII regions. 

5.5. Gas Densities in the Outer Disks: Sub-threshold Star Formation 

We examined the HI gas density at the star forming region that is furthest from the 
center for each of the five galaxies with HI maps. Using gas surface density radial profiles, 
we calculated the azimuthally-averaged HI density in Mq pc~^. These values are presented 
in Table O There is a range in gas surface density at the radius of the farthest star-forming 
region, 'LgasiRuv)) from 1.4-4.4 (1.9-5.9, corrected for He) Mq pc~^. Schaye (2004) suggests 
that the threshold for star formation is a surface density of 3-10x10^° HI atoms cm~^, or 
roughly 2.4-8 Mq pc~^. Hunter et al. (2001) find that 2-8 Mq pc~^ is required for star 
formation as defined by HII regions in NGC 2366. The values of T,gas{Ruv) found here 
are at best at the low end of these previous limits in three of the galaxies. The furthest 
star-forming regions are not incredibly massive, but high masses are not expected for single- 
member populations (i.e., furthest regions), according to the size of sample effect ( §5.6p . 
Still, there is star formation at low average gas column densities. One possibility is that 
there is a large amount of cold molecular material in outer dwarf disks. Galliano et al.(2003) 
find that in NGC 1569 there is a significant millimeter excess, perhaps indicating extensive 
cold (5-7 K) dust that may make up 40-70% of the total dust mass of the galaxy. If this dust 
is real, then this is perhaps evidence of pervasive cold molecular material. Alternatively 
and more likely, the star-forming regions may be at places where the local gas density is 
significantly higher than the average surroundings (Schaye 2004; Elmegreen & Hunter 2006). 
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The problem then is to determine the mechanism of cloud formation in sub-threshold gas. 

Most of the far outer star-forming regions in spiral galaxies shown by Thilker et al. 
(2007) are in spiral arms. Gaseous spirals can propagate much further out in a disk than 
stellar spirals because the gas spirals do not get absorbed at the outer Lindblad resonance 
like the stellar spirals do. Thus there is a natural explanation for cloud formation in outer 
spiral disks that is connected with compression from gaseous density waves. The far-outer 
regions in dwarf Irregular galaxies are not aligned in spiral arms, however, so this explanation 
does not apply here. There has to be another source of compression for Irregulars that is 
consistent with the continuation of the exponential disk. This exponential constraint would 
seem to rule out random extragalactic cloud impacts or minor mergers, which would trigger 
star formation in a more irregular way. There is also little evidence for supersonic turbulence 
in far-outer disks, so the compression from that would seem to be small. For example, the 
HI hue widths in outer disks could be largely thermal at several thousand degrees Kelvin, 
i.e., from a nearly pure warm phase (Young & Lo 1996, Young et al. 2003, de Blok & Walter 
2006). 

An interesting possibility is that star formation is triggered by gravitational instabilities 
in the outer disk, just like it is thought to be triggered in the inner disk. This can occur in sub- 
threshold gas if the angular momentum from a growing perturbation is removed. The Toomre 
Q threshold applies only when perturbations conserve angular momentum (Elmegreen 1987). 
With magnetic tension or viscosity that can remove angular momentum from a growing per- 
turbation as it tries to spin up, the disk is always unstable even at large Q. What changes 
is the growth time, which becomes long at large Q, in rough proportion to the viscous time 
or magnetic damping time (Elmegreen 1991; Gammie 1996). If the viscous and damping 
times are approximately equal to the crossing time through the thickness of the disk, which 
is true for equipartition magnetic fields and turbulent damping, then there should be no 
sharp transition at the Q — 1 radius. The instability growth time would increase steadily 
with radius in proportion to the disk thickness crossing time, and this would always give 
some level of instability for cloud formation. Previous observations of a threshold, based 
on a drop in Ha emission, would then be an artifact of low emission measure and Lyman 
continuum photon escape, as discussed above. 



5.6. Size-of-SampIe effect: Region Masses 

Large samples of star-forming regions are more likely to contain rare massive members 
than are small samples. This should be true for any subsample of a collection too. If the 
mass distribution function of star-forming regions is dN/dM oc M~^, which is usually the 
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case, then the mass of the most-massive region should increase hnearly with the number 
of regions. Previous discussions in this paper illustrated numerous examples of this size-of- 
sample effect. Figures [TS] and [12] show two results related to the size of sample effect: a 
trend in the maximum mass of a star-forming region in a radial interval versus the surface 
density of regions (Fig. [T5|l and a trend versus the total number of regions (Fig. [T6ll in that 
interval. In both cases, there is a positive correlation such that the most massive regions are 
also where the surface density and total count are highest, which tends to be the inner part 
of a galaxy. The clearest correlations have a slope of about unity. This figure illustrates an 
effect where the star-forming regions in the outer part of a galaxy tend to be lower mass 
than those in the inner part, not because of a physical change with radius, but because the 
total count of regions in the outer part is smaller. 

We are also interested in whether the size-of-sample effect applies when comparing one 
galaxy to another. Figure [T7] shows the maximum region mass plotted as a function of 
the total number of regions extrapolated down to a mass of 10'^ Mq, which is the lowest 
observable minimum mass for all of our galaxies. If the observable minimum Mjam for a 
particular galaxy is higher than this, then the extrapolated count is taken to equal the count 
higher than Mmin multiplied by the ratio Mq. This would be the extrapolated 

count in the case when the mass function is oc M~^. The figure shows a nice correlation 
with a slope of unity, as expected for the size-of-sample effect. Thus galaxies with more 
massive star-forming regions are the ones with more total regions (e.g. Billett, Hunter, & 
Elmegreen, 2002; Whitmore 2003). The galaxies have these massive regions only because the 
star- formation process is able to sample further into the tail of the region mass distribution 
function. This result implies that star formation operates with universal physical processes 
that are stochastic in nature. Of course, the specific region masses are determined by these 
processes, but the mixture of possibilities is so large that the mass distribution function is 
essentially random when it is observed as an ensemble. 

The one exception in our study seems to be the super star cluster in NGC 2366, which 
is relatively unaccompanied by lower mass clusters ( §A.6p . This giant cluster could either 



have been formed by a completely different process, one that does not give the same 
mass distribution function as the main process, or it could have grown to its current mass 
by coagulation of lower-mass clusters. Models of cluster coalescence are in Elmegreen et al. 
(2000) and Fellhauer & Kroupa (2002). 
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6. Summciry 

With the small size of dwarf galaxies comes low surface brightnesses and low densities. 
Yet despite these apparent impediments to star formation, we see dwarf galaxies actively 
forming stars, sometimes in dramatic starbursts. We also see star formation in the far outer 
parts of dwarf galaxies, where the conditions are even more extreme. Here we identified star- 
forming regions in nearby dwarf galaxies on UV images in order to study their properties as 
a function of galactocentric radius and from galaxy to galaxy. 

We measured the ages and masses of the more prominent star-forming regions in 11 
nearby dim galaxies. We were able to define 1623 useable regions on GALEX NUV images. 
We see star forming regions that run the gambit from tens to millions of solar masses, and 
we were able to find regions as young as our lower limit of several million years to our 
completeness limit of 500 Myrs. Within our sample of galaxies we see stark differences from 
one to the next. For example, our closest galaxy has at least 713 star-forming regions, while 
our next closest galaxy has only 5. Still, for all of our galaxies, the maximum region mass 
scales linearly with the total number of regions, and within each galaxy, the maximum mass in 
a radial annulus scales linearly with the local number in that annulus. Both trends indicate 
that star formation mass is independent of local environment and even the surrounding 
galaxy mass, but is given by random sampling from a universal mass function that is of the 
usual form dN/dM oc 

We find that our galaxies fall into three categories describing cluster formation activity 
over the past 500 Myrs: star formation with gaps, nearly constant star formation, and little- 
to-no star formation. The cases with gaps include some with occasional spikes in the star 
formation rate, and others with a cessation of star formation in the last 10 to 25 Myr. 

There is no consistent, measurable radial dependence for ages and masses of star- forming 
regions. The exceptions, DDO 50 and NGC 2366, show a trend towards lower masses at 
greater radii, but even there we see massive regions on the edges of the optical galaxies. 
Generally there is no indication that dwarf Irregular galaxies are less capable of producing 
massive star-forming regions in their outer disks than in their inner disks. Still, the regions 
in the outer disk often have no discernable Ha emission. This is presumably the result of the 
low density there, which makes Ra in the vicinity of massive stars very difficult to observe. 

These little galaxies have in the past presented many quandaries as to how something 
so small and sparse can form stars, and now it looks like they are forming stars in areas even 
more extreme than was previously believed possible from other surveys. This is consistent 
with a picture in which localized processes form stars in most regions of dwarf galaxies. 
At locations where the average gas surface density is much smaller (xlO) than the fiducial 
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critical value, we still see star formation. The detailed reasons for this are unknown. There 
are no spiral distributions for this star formation as there are in the outer parts of large 
spiral galaxies, so we cannot identify spiral wave compression as a primary cloud formation 
mechanism. It could be that all of the outer galaxy star formation and much of what 
happens in the inner region too, is triggered locally, by stray supernova, for example, or 
by the turbulence that supernovae and other energy sources drive. It could also be that 
large-scale gravitational instabilities operate everywhere, even in sub-threshold gas, because 
of an easy transfer of angular momentum away from growing density perturbations. In that 
case, the collapse rate could be smaller than the local dynamical rate, more like the angular 
momentum loss rate, but the final result for the scale and properties of star formation could 
be the same. High resolution observations of the gas in these regions should clarify this 
picture. 

The star forming regions in our survey also offer some insight into the evolution of 
measureable mass with age. The distribution of regions on a logM versus logt diagram 
is somewhat uniform over age t for equal intervals of mass M. This uniformity suggests 
that the regions are being removed from our view systematically with age. This loss is 
consistent with a model in which each individual region mass is decreasing as 1/t. It is also 
consistent with another model in which each region mass is constant until that region is 
suddenly disrupted, and the disruption probability is proportional to 1/t. Both possibilities 
could happen simultaneously too. We cannot say anything about the regions that have 
been disrupted, but we note from a trend in the age dependence of surface brightness that 
measured masses may in fact be decreasing as 1/t. This result is illustrated most clearly in 
a diagram of fading-corrected absolute magnitude versus the log of the age. A correlation 
in this diagram has a slope consistent with M oc 1/t. The implication is that the outer 
parts of star-forming regions are getting lost from our view as they slip below the surface 
brightness limit of the survey. This is a loss of mass over and above stellar evolutionary 
effects (supernovae and winds) and in addition to evaporation. We are investigating all of 
these mass loss processes in more detail in a second paper. 
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Facilities: GALEX, Lowell Observatory 

A. Individual galaxies 

The UV-bright knots identified in each galaxy are outlined in Figure [181 Here we discuss 
the properties of the ensemble of regions in each galaxy. 

A.l. NGC 6822 

NGC 6822 is the closest galaxy in our survey (0.5 Mpc). Because of this, we were able 
to identify 800 individual star-forming regions, 713 of which made it into our final data set. 
With this galaxy we were able to trace the star forming regions to much lower luminosity 
levels than in any other galaxy. A problem arises because of the galaxy's large angular size 
(about 15'xl5'), which led to field star contamination at this low galactic latitude. 

NGC 6822 has star-forming regions ranging in mass from 15 Mq to 3 x 10^ Mq, and 
ranging in age from several Myrs to well beyond the 500 Myrs that is emphasized in the UV. 
The lower mass limit of a star-forming region that should still be visible after 500 Myrs is 
5 X 1O^M0 in this galaxy. Figure H] suggests that star formation has occurred in two or three 
major episodes. These episodes appear as peaks in the number of all clusters detected and 
in the number of clusters detected with masses larger than our adopted lower mass limit. 
They also appear as peaks in the mass of the most massive cluster as a function of age. The 
episodes occur for logarithmic ages between the present and 7.5, between 8.12 and 8.78, and 
after 9.28. We refer to these as the 7 Myr, 100 Myr, and 1 Gyr bursts. 

In the most recent star formation event at 7 Myr, there are 3 clusters more massive than 
5 X 10^ Mq. The most massive cluster in this event is 2 x 10^ Mq. In the event at 100 Myr, 
there are 7 regions of mass > 5 x 10^ Mq, with the most massive of these having a mass of 
7 x IO^Mq. In the older burst at IGyr, which may not be complete in time and is certainly 
not complete to our adopted lower mass limit, we see 54 regions of mass > 5 x IO^Mq, with 
the most massive having a mass of 3 x IO^Mq. The oldest age burst is old enough that all 
of the clusters below 5 x IO^Mq have faded, and we see this in the lower mass limit of the 
regions we identified. 

We can estimate the apparent relative formation rates of star-forming regions from the 
counts of regions on the time axis. For the 7 and 100 Myr bursts, the apparent region 
formation rates with M > 5 x 10^ Mq, obtained from the counts per unit time interval, 
are 0.1 and 0.01 massive regions per Myr, respectively. These are the apparent and not the 
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real formation rates because many clusters have been lost from their original mass intervals 
through disruption and evaporation. 

For a mass function of the form dN/dM oc M~^, the total number of initial regions 
should be equal to the ratio of the maximum region mass to the minimum region mass. 
This is one aspect of the size-of-sample effect. If we consider a minimum detectable mass to 
be 5 X 10^ Mq, then the number of regions above this mass should be proportional to the 
maximum region mass if the disruption probability is independent of mass (e.g., as observed 
by Chandar et al. 2009). From the pure counts of clusters above 5 x 10^ Mq, the proportion 
in the two most recent bursts is 1:2.3, and from the maximum masses, the proportion is 
1:3.5. These proportions are close enough to each other to agree with the predictions of the 
size of sample effect, considering the small numbers of massive clusters, (fractional errors in 
the counts are 60% and 38%, respectively, from the inverse square roots of the number of 
massive regions). 

Fall et al. (2005) consider the mass-age diagram for the Antenna galaxy and note that 
generally the number of clusters above a certain mass in equal log intervals of time is about 
constant with time, which means that the number of clusters remaining per unit age interval 
decreases as 1/age. In our sequence of ages, the counts per unit time, 0.1 and 0.01 regions 
Myr~^ are indeed roughly proportional to the inverse of their ages (7 and 100 Myr). This 
implies that clusters of all masses are becoming increasing lost over time. Such loss could be 
from a combination of the two processes as discussed in §5.21 and from the surface brightness 
effects discussed in §5.11 

From logt = 7.9 to logt = 8.1 and from logt = 8.8 to logt = 9.3, periods of 46 Myr 
and 1300 Myr, there are lulls in the star formation rate. In the first interval, no regions of 
M > 1 X 10^ Mq form and the maximum mass drops down to 1 x 10^ Mq. In the second 
interval, no regions were formed as massive as the periods on either side, but the maximum 
mass is 1 X 10^ Mq. In the bursts before and after these lulls, the maximum cluster masses 
were higher by factors of ~ 7 — 70 and the total numbers of clusters remaining were higher 
by factors of ~ 1 — 10. Thus we are seeing oscillations in the cluster formation rate by a 
factor of order ~ 10 in both the cluster counts and the maximum cluster mass. This type of 
variation is consistent with Lee et al.'s (2007) prediction. 

The time intervals of the peaks and lulls decrease toward the present time. This appears 
to be the result of star formation with a clustered or hierarchical structure in time (for models 
of such star formation, see Elmegreen & Scalo 2006). Cluster bursts of short duration that 
occurred long ago would appear blended in time because our time resolution cannot resolve 
the small subbursts, while cluster bursts of the same intrinsic width that occurred recently 
can still be resolved. The time resolution increases in proportion to a power of the cluster 
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age because the cluster colors evolve with a power law time dependence. 

Figures [3] and H] suggest that the recent formation of massive regions in NGC 6822 is 
more prominent in the outer part of the disk, beyond Rho- The left-hand part of Figure H] 
also suggests that most of the young regions, regardless of mass, are outside Rna- In fact 
468 of this galaxy's 713 star-forming regions fall outside of Rho, which is far more than we 
see in any other galaxy in our survey. This is consistent with the distribution of blue stars 
(Komiyama et al. 2003). 

A. 2. IC 1613 

IC 1613 is our third nearest galaxy at 0.7 Mpc, and second most populous galaxy with 
342 useable star-forming regions. Because it is nearby, we were able to trace the star-forming 
regions to masses as low as ~ 30 Mq. The minimum mass for which we should see all 
currently definable star- forming regions from the last 500 Myrs is approximately 3 x IO^Mq, 
as indicated by the dashed horizontal line in Figure [31 The maximum mass formed is of 
order 3 x lO^M©. In IC 1613 there has been more-or-less continuous star formation over the 
last 1 Gyr. The apparent formation rates of massive star-forming regions (> 3 x 10^ Mq for 
IC 1613) from log age of 7.0 to 8.5 is 0.03 regions Myr~^. 

IC 1613 has an Ha profile that ends before the UV star-forming regions end. The Ha 
only extends to a radius of 7.5' from the center of the galaxy, and we see star-forming regions 
in the UV out to a radius of 9.3'. There are far more massive regions inside Rna than outside. 
The clusters beyond Rna tend to be lower mass with a few notable exceptions among the 
oldest identified clusters. 



A.3. WLM 

Wolf-Lundmark-Melotte (WLM) was the third most populous galaxy in our survey with 
191 regions marked, 165 of which made it into our final data set. Its distance was typical at 
1.0 Mpc. This galaxy is distinct in having the oldest average star-forming regions in all our 
sample. The cluster formation rate appears to have been relatively constant until it stopped 
about 22 Myrs ago. The minimum detectable mass at 500 Myrs is about 1 x 10^ Mq. 

This galaxy has a large fraction of star- forming regions beyond Rua'- 31 out of the total 
165. UV regions are seen to a radius of 7.4', while Ha drops out beyond 4.2'. In Figure [31 
there appears to be a slight drop in the mass beyond Rho, amounting to a factor of ~ 3. 
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A.4. DDO 50 

DDO 50 was the second most distant galaxy in our survey at 3.4 Mpc, but we were 
still able to define 169 regions, 139 of which made it into the final data set. Had this 
galaxy been located as close as WLM, from the number of massive regions we see now and 
our experiments with WLM images, we might expect to have identified 520 regions instead 
of only 169. Cluster formation in this galaxy appears to have been relatively constant 
until about 9 Myrs ago. However, there is a small enhancement in the numbers of massive 
clusters 23-30 Myrs compared to the period of time before that. The most massive cluster 
we identified is 1 x 10^ Mq . 

Because Ha extends nearly as far as star formation seen in the UV (4.1' for Ha and 
5.1' for the UV), we found only 10 regions beyond Rua- The outer regions sample the full 
distribution of masses despite the small sample size, but they tend to be a little lower in 
mass on average compared to those in the rest of the galaxy. 

A. 5. DDO 75 

DDO 75, also known as Sextans A, is another populous galaxy (119 useable regions) 
at a moderate distance (1.3 Mpc). One of the unusual characteristics of this galaxy is 
that its azimuthally-averaged \^-band and 4. 5- micron surface brightness profiles do not fall 
off at a constant exponential rate. Instead, we see a flat profile in the inner 1.8' radius 
(0.7 kpc), followed by a normal exponential decay curve beyond that (Hunter & Elmegreen 
2006). To investigate whether the break in the profile represents a change in star formation 
characteristics, we compare the regions interior and exterior to the break in Figure [61 We 
plot mass versus age for the inner and outer disks in the top part of the figure, the number of 
regions in radial intervals of 0.25 arcmin versus radius in the bottom left, and the average My 
versus radius in the bottom right. In none of these plots do we see a significant change in the 
star- formation characteristics at the surface brightness break. There is a peak in the surface 
density of star-forming regions at the break, but no obvious change in region properties. 
This result implies that the processes creating the break in the exponential \^-band surface 
brightness profile do not affect the masses or the ages of the star-forming regions. Most of 
the regions, and all of the massive regions, formed inside -Rho, so in Figure [3] we present a 
comparison of the inner half of the star-forming regions to the outer half. 

DDO 75 appears to have had a fairly constant level of star formation for the last several 
hundred million years until about 9 Myrs ago. In the past 9 Myrs only one cluster has 
been formed. Furthermore, clusters older than about 200 Myrs are missing except for a few 
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especially massive ones: clusters up to 4 x 10^ M© in mass. 

The apparent formation rate above the detectable mass limit of 2 x IO^Mq for ages 
between 9 Myr and 200 Myr was^ 0.08 Myr~^. The maximum cluster mass is 2 x 10^ Mq. 
This steady low rate is consistent with the global view of the galaxy from iJ5'T images and 
color-magnitude diagram fitting (Dolphin et al. 2003). 

A.6. NGC 2366 

NGC 2366 is among the more distant galaxies in our survey, at 3.2 Mpc, but despite 
this we were able to define 66 regions, 58 of which were useable. This galaxy also has JHK 
data. NGC 2366 contains the most massive star cluster in our combined sample, one with 
M ~ 7 X 10^ Mq. a second cluster has a mass of 2 x lO^M©. These two regions are I and 
II, respectively, in Drisscn et al. (2000). Drissen et al. split cluster I into two subclusters, A 
and B, but this divide was not discernible in our GALEX images, most likely because lA is 
still completely enveloped in dust and gas. Drissen et al. find ages of 10 Myrs for cluster II, 
2.5-5 Myr for cluster IB and an age less than 1 Myr for cluster lA. Our ages arc significantly 
older for both regions, 40 Myrs for cluster II and 700 Myr for cluster I. Our older ages may 
be due in part to our underestimating the reddening, which leads to intrinsic colors that are 
too red. Another factor is the relative sizes of the apertures we used. Our larger apertures 
encompassed the entire complexes, including outlying cluster members. Billett et al. (2002) 
also looked at cluster I, or at least two smaller chunks of it. They re-examined IB from 
Drissen et al. and what they called cluster 1, both of which were given radii of 4.6 pc - much 
smaller than the 200 pc radius used for the whole complex here. They also find My = —9.5 
for cluster IB and My = —8.5 for cluster 1, as compared to My = —14.4 for the whole 
complex here. 

The minimum mass for which we expect completeness up to 500 Myrs of age is 1 x 10"^ Mq, 
about the same as for DDO 50. The number of clusters more massive than this is also about 
the same in these galaxies (21 for NGC 2366 and 31 for DDO 50), and the most massive 
cluster is extremely massive in both: 7 x 10^ Mq in NGC 2366 and 1 x 10^ Mq in DDO 50. 
A surprising difference is that NGC 2366 has far fewer low mass clusters in the 10 Myr to 
100 Myr age interval, 19 smaller than 10^ Mq in NGC 2366 compared to 93 of this mass m 
DDO 50. The cluster mass function is apparently flatter in NGC 2366, or more concentrated 
in its most massive two clusters. This difference could be the result of cluster coalescence in 
NGC 2366, which would explain both the lack of low mass clusters and the excessive mass 
of the largest cluster in that galaxy. Cluster coalescence was shown to be reasonable for 
another 10^ Mq cluster in NGC 6946 (Elmegreen, Efremov & Larsen 2000). 
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This galaxy has only one NUV region beyond Rho- We present instead a comparison 
of the inner half of the star-forming regions to the outer half in Figure [31 We see several 
low mass, young star-forming regions in the outer half that are absent in the inner galaxy, 
and in the very outer galaxy (past 4 kpc radius) we see no massive regions. NGC 2366 is a 
good example of a galaxy with a gradient in the star formation region mass. Inside 4 kpc, 
the > IO^Mq regions are uniformly distributed. 

A. 7. DDO 70 

DDO 70 is a small galaxy (47 useable regions) at a moderate distance (1.3 Mpc). It has 
had steady, slow star formation over the last 200 Myrs, with peak masses of only ~ 8 x IO^Mq. 
Like DDO 75, DDO 70 has formed only a few clusters over the past 10 Myrs and only a few 
clusters before 200 Myrs ago. There are no star-forming regions found beyond Rho, and so 
no comparison can be made between the galaxy interior and exterior to Rna- A plot of age 
versus mass is presented in Figure [3] with regions color-coded according to the inner 50% 
and the outer 50% of regions. 

A.8. DDO 216 

DDO 216 is a very faint little galaxy. Despite its relatively low distance (0.9 Mpc), we 
were able to identify only 34 regions, 25 of which have made it into the final data set. This 
galaxy has a double exponential y-band surface brightness profile: The profile interior to a 
radius of 5.4' drops at one rate, and the profile exterior to that radius falls off at a faster 
rate. However, there are no UV-identified star-forming regions beyond a radius of 5.4'; our 
furthest region was only 3.0' radius from the center of the galaxy (deprojected distance). 

DDO 216 has had very little star formation over the last 20 Myrs, but before that there 
was a steady stream of modest star formation back to at least 500 Myr. The lack of young 
star clusters accounts for the small cutoff radius, Rho (which is 1.3'), and consequently most 
of the UV-detected star-forming regions are exterior to Rho- 14 of 25 are beyond Rua- The 
only difference in the regions inside and outside Rho is that those outside extend to older 
ages. We see no difference in the minimum age or maximum mass of clusters inside and 
outside of Rhu, as shown in Figure [31 
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A.9. DDO 210 

DDO 210 is a tiny little wisp of a galaxy. It is the second faintest galaxy in our 
sample at an My of —10.9, significantly fainter than some of our super massive star clusters. 
Furthermore, this galaxy has no HII regions. Because of its extreme faintness, and despite 
its relative proximity (0.9 Mpc), we were able to find only 24 regions in the NUV image, 
and only 9 of those made it into the final data set. Of these nine regions, one is less than 10 
Myrs old, and 6 are 100 Myrs old or older. The lack of Ha for the youngest regions reinforces 
our suggestion elsewhere in this paper that Ha is too faint to see despite the likely presence 
of ionized gas. Figure [3] separates the regions into the inner 50% and outer 50%. 

A. 10. LGS3 

LGS 3 has only 7 definable regions, 5 of which were useable. This was the second closest 
galaxy in our survey at a distance of 0.6 Mpc, so we should be able to see any star formation 
that might have been there. Despite these favorable conditions we still see nearly no star 
formation over the last billion years. This galaxy has no measurable HII regions. 

A.ll. MSldwA 

For this tiny little galaxy we were able to find only three star-forming regions, two of 
which were useable. MSldwA is the most distant galaxy in our survey at 3.6 Mpc. The two 
regions we do see are reasonably large, but between the two of them, they compose almost 
the entire galaxy in the NUV. This is another galaxy without any measurable HII regions. 
As we have only two data points, statistical comparisons are not possible. 
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Fig. 1. — Comparison of brightness distributions at different distances. Tlie number of 
regions as a function of My of the regions is shown for DDO 50, which is at a distance of 
3.4 Mpc, and for NGC 6822, which is much closer at 0.5 Mpc. The affects of distance and 
blending on our detection limits is evident in the larger percentage of fainter regions detected 
in NGC 6822. 
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Fig. 2. — Plots of mass vs. age, mass vs. galactocentric radius, and age vs. radius for WLM. 
The nominal distance of WLM is 1 Mpc, and the panels labeled as 1 Mpc refer to regions 
identified and measured on the original images. The panels labeled as 3.4 Mpc refer to 
regions identified and measured on images that have been smoothed to the resolution we 
would have if the galaxy were located at 3.4 Mpc. Regions outside the extent of Ha are in 
red. The uncertainties in the ages and masses have been left off in order to make the plot less 
crowded. At 3.4 Mpc, we see fewer regions with more mass and shallower radial gradients, 
as expected, but general trends remain. 
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Fig. 3. — Plots of mass vs. age, mass vs. galactocentric radius, and age vs. radius for the 
galaxies in our survey. Regions outside the Ha extent or in the outer half of all regions are in 
red, as indicated by the labels. The horizontal dashed lines in the Mass- Age plots are mass 
limits for completeness to an age of 500 Myrs. The slanted dashed lines are fits by eye to 
the upper envelopes of the cluster distributions. The slanted solid line shows the slope for a 
fading relationship in which the minimum observable mass scales as log Mass oc 0.69 log Age. 
The vertical dashed lines correspond to radii that separate the inner and outer parts, as 
defined by the colored symbols. 
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Figure [3] (continued) 
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Figure [3] (continued) 
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Figure [3] (continued) 
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Fig. 4. — Total number of regions in our survey and the number larger than the minimum 
observable mass (as indicated by each label) shown as a function of region age with an age 
interval of A log Age = 0.2. On the right, the maximum region mass in each age interval is 
shown. Red histograms are for outer galaxy star formation, in the sense defined by Figure 

El 
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Figure H] (continued) 
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Fig. 5. — Mass (left) and age (right) distributions for the inner and outer halves of star- 
forming regions in NGC 6822 {top two panels of figure), DDO 50 {middle two panels), and 
NGC 2366 {bottom two panels). The inner half of the regions is in the top panel for each 
galaxy, and the outer half is in the bottom panel. 
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Fig. 6. — Distributions of star formation interior and exterior to the surface brightness break 
Rsr in DDO 75. In the upper left are plotted the masses and ages of regions inside the break, 
and in the upper right are plotted the masses and ages of regions outside the break. The 
lower left plot is the number of star-forming regions per 0.25', and the lower right plot is the 
average My for the star-forming regions in 0.25' bins. The vertical line in the lower panels 
marks Rbt- 
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Fig. 7. — Average NUV surface brightness and number density of all the identified regions 
shown as functions of galactocentric radius for each galaxy with a large number of regions. 
The vertical dashed lines correspond to radii that separate the inner and outer parts, as 
defined in Figure [3l 
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Figure [7] (continued) 
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Figure [7] (continued) 
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Fig. 8. — Filling factor of star-forming regions shown versus galactocentric radius. The 
filling factor equals the sum of the region areas divided by the area of the radial annulus. 
The trends here follow the radial trends in luminosity density and counts because the UV 
region sizes are about constant with radius. The vertical line separates the inner and outer 
regions as defined in Figure [3l 
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Fig. 9. — Radius of each region plotted versus the galactocentric radius. The region sizes 
are estimates based on the polygons used to visually identify each star-forming region in the 
NUV GALEX images. They are not accurate measures of region size. Still, they should 
reveal systematic differences of a factor of two or more in the sizes for the inner and outer 
parts of the galaxy (as delineated by the vertical lines), and no such differences are observed. 
The regions in DDO 50, NGC 2366, and MSldwA are significantly larger compared to regions 
in other galaxies in our sample. This is likely due at least in part to the larger distances 
of these galaxies, and the difficulty in defining star-forming regions, especially faint, small 
regions, to the same degree of precision as was possible in closer objects. 
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Figure [9] (continued) 
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Fig. 10. — NUV surface brightness of each star-forming region in four galaxies, versus the 
distance from the galactic center. The surface brightness decreases with distance for DDO 
50 partly because of a selection effect where regions have to be brighter than the surrounding 
disk in order to be picked for the survey. This effect is also slightly present for IC 1613 and 
WLM, but not obviously for NGC 6822. 
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Fig. 11. — NUV surface brightness of each star- forming region in four galaxies, versus the 
log of the region age. Dashed lines indicate the fading trend from the evolution of a single 
stellar population model. Each galaxy has a lower limit to surface brightness that depends 
on exposure time but not age. The top left panel suggests a loss of surface brightness purely 
from stellar evolution fading, while the other panels suggest some additional loss because 
the distributions are steeper than the fading trend. 



-54- 



140 
120 
100 
80 
60 
40 



o 20 



^ 180 
5 



1 ' , 1 ' 1 

7NGC 6822 " 

• 

• • 

. • • . . 

': \ ^ \ ^ \ ': 


1 ■ 1 ' 1 
~IC 1613 ~ 

• 

• 

• • 

': H-H \ ^ \ ; 


• 

• 

I.I.I 


rODO 50 . ~ 
i • • 

•* • 

— • — 

1 . 1 . 1 ■ 



9 7 
log Age (yrs) 



Fig. 12. — NUV region radius is shown versus log age. The regions have about the same 
size for all ages, indicating no expansion of their visible outer boundaries. 
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Fig. 13. — The fading-corrected absolute magnitude is plotted versus the log of the region 
age. The dashed line has a slope of +2.5 and traces the outline of dense points, which show 
an upward trend. This trend suggests each region has a mass that varies as 1/age. 
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Fig. 14. — A schematic model showing projected profiles (Gaussians) through two clusters 
with the same total mass but different masses above the surface brightness limit. The green 
rising curves are the cumulative masses above the surface brightness limit, integrated from 
left to right over the profiles. The mass above the surface brightness limit is smaller for 
the fainter region, even though the total mass is the same. This model illustrates how 
observable regions can appear to lose mass over and above any effects from stellar evolution, 
evaporation, and disruption. 
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Fig. 15. — Maximum cluster mass in equal intervals of galactocentric radius plotted versus 
the number of regions per unit area in those intervals. Bins are in 200 pc increments of 
radius. 
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Fig. 16. — Maximum cluster mass in equal intervals of galactocentric radius plotted versus 
the total number of regions in those intervals. Bins are in 200 pc increments. The outer 
parts of most galaxies tend to have the smallest numbers of regions. There is a trend from 
the size of sample effect with a slope of unity on this diagram. 
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Fig. 17. — Maximum region mass in each galaxy, regardless of position or age, plotted versus 
the total number of regions extrapolated to a mass of 10^ Mq. The power law correlation 
with slope 1 follows from the sizc-of-sample effect in a population of star-forming regions that 
have a mass function dN/dM oc M~^. This correlation suggests that larger and more active 
dwarf galaxies have more massive clusters entirely because of sampling statistics, without 
any physical differences in the detailed star formation mechanisms. 
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Fig. 18. — NUV images of each galaxy with individual star-forming regions outlined. For 
regions outlined in green, the background was determined from a surrounding annulus. For 
regions outlined in red, the background was determined by sampling areas around the star- 
forming region that seemed to represent the underlying stellar disk. 
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Figure [18] (continued) 
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Figure [18] (continued) 



Table 1. Galaxy Sample^ 









D 






Rd'= 


P.A.d 


Incl.° 


Galaxy 




Other Names 


(Mpc) 


E(B - V)f^ 


Mv 


(arcmin) 


(dcg) 


(dcg) 


DDO 50 . 


PGC 


23324,UGC 4305,Holmberg II,VIIZw223 


3.4 


0.02 


-16.6 


1.11 


18 


46 


DDO 70 . 


PGC 


28913,UGC 5373,Sextans B 


1.3 


0.01 


-14.1 


1.26 


88 


57 


DDO 75 . 


PGC 


29653,UGCA 205,Sextans A 


1.3 


0.02 


-13.9 




41 


33 


DDO 210 


PGC 


065367,Aquarius dwarf 


0.9 


0.04 


-10.9 


0.63 


-85 


66 


DDO 216 


PGC 


71538,UGC 12613,Peg DIG 


0.9 


0.02 


-13.3 


1.68 


-58 


69 


IC 1613 . 


PGC 


3844,UGC 668,DDO 8 


0.7 


0.00 


-14.6 


2.87 


71 


37 


LGS3 .... 


PGC 


3792, Pisces dwarf 


0.6 


0.04 


-9.4 


1.13 


-3 


64 


MSldwA 


PGC 


23521 


3.6 


0.02 


-11.7 




86 


45 


NGC 2366 


PGC 


21102,UGC 3851,DDO 42 


3.2 


0.04 


-16.7 


1.37 


32 


72 


NGC 6822 


PGC 


63616,10 4895,DDO 209,Barnard's Galaxy 


0.5 


0.21 


-15.2 


3.92 


24 


40 


WLM ... 


PGC 


143,UGCA 444,DDO 221 


1.0 


0.02 


-14.4 


1.97 


-2 


70 



^Quantities taken from Hunter & Elmegreen (2006). 
''Foreground reddening from Burstcin & Heilcs (1984). 
'^Disk scale length determined from V-band images. 

''Galactic position angle of the major axis determined from V-band images. 
''Inclination of galactic disk. 
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Table 2. GALEX Observations 





FUV Exp (s) 


NUV Exp (s) 


Tile Name 


DDO 50 . 


1674 


1521 


NGA_HolmbergII 


DDO 70 . 


1076 


966 


NGA_SextansB 


DDO 75 . 


1663 


1512 


NGA_ScxtansA 


DDO 210 


1669 


1472 


GIl_047107_DDO210 


DDO 216 


1557 


2665 


MISDR2_28664_0746 


IC 1613 . 


1676 


1494 


NGA_IC16134 


LGS3 .... 


1673 


1487 


NGA_LGS3 


MSldwA 


1674 


1521 


NGA_HolmbergII 


NGC 2366 


2889 


2617 


NGA_NGC2366 


NGC 6822 


4483 


5132 


NGA_NGC6822 


WLM ... 


1423 


1294 


NGA_WLM 



Table 3. Region Photometry 



R.A. Dccl. Trcg^ log M'* logAge'^ 

{J2000) {J2000) (kpc) RIRd°' (pc) Al"lV!71/° O" (FUV - NUV)o u Mv O" (B - V)o cr {U - S)o u (Mq ) (yrs) RMS'^ 



DDO 50 



1 


8 19 09.6 


70 45 19 


1.96 


1.78 


108.1 


-9.9 


0.01 


-0.08 


0.03 


-9.5 


0.0 


-0.01 


0.02 


-0.94 


0.03 


4.03 


7.35 


0.38 


2 


8 18 49.7 


70 44 60 


2.99 


2.72 


81.4 


-8.6 


0.02 


-0.44 


0.04 


-8.5 


0.0 


-0.02 


0.04 


-1.00 


0.05 


3.55 


7.08 


0.38 


3 


8 19 12.1 


70 43 07 


0.55 


0.50 


191.6 


-12.1 


0.00 


-0.17 


0.01 


-12.2 


0.0 


0.21 


0.00 


-1.14 


0.00 


5.08 


7.27 


0.55 


4 


8 19 27.8 


70 42 21 


2.65 


2.41 


139.4 


-11.3 


0.01 


-0.22 


0.01 


-11.0 


0.0 


0.05 


0.01 


-1.20 


0.01 


4.53 


7.03 


0.55 


5 


8 19 28.6 


70 42 30 


2.65 


2.40 


64.4 


-8.7 


0.02 


-0.15 


0.04 


-8.0 


0.0 


-0.05 


0.04 


-1.09 


0.05 


3.41 


7.20 


0.49 



Note. — Table|3]is published in its entirety in the electronic edition of the Astronomical Journal. A portion is shown here for guidance regarding its form and content. 

'^Distance of the star-forming region from the center of the galaxy in the plane of the galaxy, using the galactic parameters given in Table [T] Rn is the disk scale length from 
Hunter & Elmegreen (2006). 

''Radius of the star-forming region. The identifying regions are polygons, and the radius is the square-root of the area of the region divided by tt. 
'^Absolute AB magnitude in the NUV filter. 

"^Region stellar mass and age from the best fit model to the colors and luminosity, and RMS. Here RMS is the rms deviation between observed and modeled colors. ^ 
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Table 4. NGC 2366 Region JHK Photometry 
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066 
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Tablc 4 — Continued 



Region 


i\/r 


(J 


[J - U)o 


(J 


TJ 

11 - 


- J\ 


a 


32 


-9.510 


1.560 


-0.031 


6.561 










33 


-10.935 


0.407 


0.534 


1.137 


0, 


.151 


2. 


,854 


34 


-9.537 


0.989 


0.485 


2.862 


0, 


.077 


7. 


,761 


35 


-9.494 


1.458 


-0.258 


7.256 










36 


-9.962 


0.887 


-0.364 


4.760 


-0, 


.062 


15. 


,992 


37 


-11.469 


0.300 


0.298 


0.992 


0, 


.119 


2. 


,648 


38 


-10.357 


0.438 


0.651 


1.129 


0, 


.146 


2. 


,794 


39 


-8.240 


2.872 


1.397 


4.743 


-0, 


.634 


20. 


,645 


40 


















41 


-11.936 


0.180 


0.617 


0.476 


0, 


.155 


1. 


,175 


42 


















43 


-8.874 


2.827 


-1.226 


28.369 










44 


-8.561 


3.035 


-0.636 


1 •■• 


1, 


.616 


20. 


,540 


45 


-7.926 


4.552 


-1.098 


41.742 


2, 


.620 


36. 


,823 


46 


















47 


















48 


-7.859 


3.774 


0.508 


10.698 


0, 


.058 


29, 


.730 


49 


-9.369 


1.381 


1.259 


2.449 


0, 


.571 


3. 


,609 


50 


-7.974 


4.131 


0.808 


9.603 










51 


-8.102 


2.600 


0.512 


7.373 


0, 


.815 


10. 


,638 


52 


















53 


-9.785 


1.041 


0.100 


3.976 










54 


-11.279 


0.329 


0.652 


0.847 


0, 


.185 


2. 


,020 


55 


-13.037 


0.089 


0.638 


0.230 


0, 


.185 


0. 


,550 


56 


-13.763 


0.046 


0.282 


0.156 


0, 


.446 


0. 


,312 


57 


-10.459 


0.419 


0.614 


1.106 


0, 


.172 


2. 


,690 


58 
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Table 5. Galactic Parameters 







n b 
IlBr 




ituv 


^half 






T.Hi{Ruvr 


Galaxy 


^^ \ 

(kpc) 


^^ \ 

(kpc) 


AT C 


(kpc) 


(kpc) 


70 past Rna 


% past UBr 


[Mq pc ") 


DDO 50 . 


4.02 




139 


5.07 


2.70 


7 




3.6 


DDO 70 . 


1.23 




46 


1.35 


0.60 


2 






DDO 75 . 


1.17 


0.67 


119 


1.30 


0.69 


2 


52 


3.4 


DDO 210 


0.00 




9 


0.20 


0.13 






4.4 


DDO 216 


0.34 


1.42 


25 


0.79 


0.36 


56 







IC 1613 . 


1.53 




342 


1.89 


0.91 


5 




1.5 


LGS3 .... 


0.00 




5 


0.25 


0.16 








MSldwA 


0.00 


0.42 


2 


0.47 


0.37 




50 




NGC 2366 


5.25 




58 


5.47 


2.79 


2 




2.2 


NGC 6822 


1.65 




713 


5.88 


2.16 


66 






WLM ... 


1.24 




165 


2.16 


0.78 


19 




1.4 



'^Radius beyond which Ha emission is no longer detected (Hunter & Elmegreen 2004). 
Galaxies with Rna of zero have no detected Hct emission. 

Radius where the F-band surface photometry changes slope (Hunter & Elmegreen 2006). 

^Number of retained regions identified on the NUV image. 

'^Radius where last star-forming region is detected in the NUV image. 

^Radius interior to which is found half of the star-forming regions identified on the NUV 
image. 

^Percentage of identified UV regions found beyond the radius Rna where Ha emission is 
no longer detected and beyond the break radius Rbv for those galaxies with breaks in their 
broad-band surface brightness profiles. 

^HI surface density at Ruv- 



